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Synopsis 

The density matrix renormalization group (DMRG) method is apphed 
to the interaction round a face (IRF) model. When the transfer matrix is 
asymmetric, singular-value decomposition of the density matrix is required. 
A trial numerical calculation is performed on the square lattice Ising model, 
which is a special case of the IRF model. 
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The renormalization group is one of the basic concepts in statistical 
mechanics. Its real-space expression — the real-space renormalization 
group^) — has been applied to various statistical models. Recently White 
proposed a precise numerical renormalization algorithm, which is known as 
the density matrix renormalization group (DMRG) method. ^'^^ This method 
has been widely used for one-dimensional (ID) quantum lattice models^"^) 
because of its numerical accuracy and portability. 

The DMRG method was originally designed for ID quantum lattice 
models. Since cZ-dimensional quantum systems are closely related to clas- 
sical systems in d+1 dimensions,®^ the DMRG method is applicable to two- 
dimensional (2D) classical systems. In this paper we apply the DMRG 
method to the interaction round a face (IRF) model,^^ which contains various 
2D statistical models such as the Ising model and the vertex models. 

The IRF model is defined by a Boltzmann weight W on each face, which 
is surrounded by four n-state spins a. The transfer matrix of the IRF model 
is given by 

2N-1 
i=l 

where and a[ are n-state spins on adjacent rows with width 2N. (Fig. 1(a)) 
We assume open boundary conditions. In general, the transfer matrix is 
not symmetric. If the Boltzmann weights satisfy a constraint — the Yang- 
Baxter relation^'^^) — the model is analytically solvable. Here, we consider 
a numerical analysis of T^"^^^ in view of the existence of unsolvable cases. 

The DMRG method for the IRF model is expressed as a renormalization 
of the transfer matrix T^'^^\ Figures 1(a) and 1(b) show the method of 
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renormalization. First of all, we assume that T^^^^ can be renormalized into 
a product form 

where Tj^^^ and T^'^ represent renormalized transfer matrices for the left 
and the right half of the system, respectively. (Fig. 1(a)) The block-spin 
variables and correspond to groups of the n-state spins {ai . . 
and {cTf^^i . . . C2jv}, respectively. We assume that and take m different 
states, where m is much smaller than n^~^. If there is a mapping from 
f{2M) f(2M+2) fQj, arbitrary M, then f ^^iv) -^^ ^^ (2) is obtained through 

successive mapping: 

rp{A) ^ 2^(4) _^ f{6) _ _ _ f{2N-2) _^ f{2N) _ 

In each step, block-spin transformations {^lCTj^} — > ^f^^^ and {cr^^^} — > 
are required. (Fig. 1(b)) The block-spin transformations are given in the 
following. Since the argument is common to both and we discuss the 
renormalization for only. 

The matrix dimension of T'(2-^) is (nm)^ for arbitrary M. If m is small 
enough, we can numerically solve the eigenvalue problem 

E T('^H^yyR^R\^L^L^dR)mL^L^RU) = ^ i^yyR^ , 

(4) 

and obtain the 'right' eigenvector ^ that corresponds to the largest eigenvalue 
In the same way, we obtain the 'left' eigenvector $ that satisfies 
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^j.(2M) ^ ^xC^M)^ j^Q^g ^Yiat if f(2^) is asymmetric, $ is not the complex 
conjugate of \E'. The density matrix'^'^-' is then expressed as a partial product 
between ^ and 

Apart from the density matrix for ID quantum systems, pl is not always 
Hermitian. Therefore we must perform the singular-value decomposition of 

PL 

E Oii\eL<T'L)PLieM^L<rL)Qi^L<rL\j) = S^jUJj, (6) 

and obtain O and Q that satisfy ^ 0{i\^LO'L)Q{^LO-]^\j) = 5^^-. We 
assume the decreasing order of the singular values: u>i > > . . .u)^^. 
The m by nm matrix 0{^'^^\^]^ai) together with the the nm by m one 
Q{^l^l\^T'^) represent the block-spin transformation from ^j^ctj^ to a new m- 
state block-spin The transformation naturally gives a relation between 

T^^^ and which is a linear transformation 

fi^^'\er<y'L\ir<yL) = 

We also obtain a similar relation between and in the same way. 

Now we obtain the renormalization processes in eq.(3). 

As one repeats the mapping in eq.(3) with the aid of eqs.(4)-(7), the 
eigenvectors of the renormalized transfer matrix tC^^) approach their fixed- 
point values. After $ and \E' have converged, we can obtain local thermo- 
dynamic quantities in the large- A/" limit. For example, the nearest-neighbor 
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spin correlation is expressed as 

when the inner product ($, is unity. In the same way, we can obtain the 
internal energy, susceptibilities, and specific heat. Observations of correlation 
functions and surface tensions are possible with the use of the 'finite chain 
algorithm' of the DMRG method.^) 

We examine efficiency of the DMRG method. As a reference, we calcu- 
late specific heat of the nearest-neighbor Ising model, where the Boltzmann 
weight is 

= exp |y ((J-CT-+1 + (j-(Ji + cT-+i(Ji+i + (Ji(Ji+i)| (9) 

for (J — ±1. In this case, the renormalized transfer matrix T*^^^-* is real- 
symmetric for arbitrary AT, and therefore the eigenvectors ^» and in eq.(5) 
are the same. The specific heat Cv is obtained from a numerical differential 
of the nearest-neighbor spin correlation (cTjCr^^.]^) in eq.(8). The numerical 
error e in Cv is obtained through a comparison between calculated data and 
the exact value. 

Figure 2 shows the calculated specific heat. The plotted data are ob- 
tained within 1024 iterations (2A^ = 2048) when m — 60. The numerical 
error e in Cv is non- negligible near Tc, partly because the numerical conver- 
gence of the bond energy {a^a^j^-^) with respect to N and m becomes worse as 
the parameter K = J/ksT approaches its critical value K* = J/ksTc. The 
calculated energy per site at T = Tg is —1.41398 when m = 60, which is com- 
parable in accuracy to a recent numerical estimate^^^ —1.419(1) by a micro- 
canonical Monte Carlo simulation; the exact value^^^ is — -\/2 = —1.414213. 

5 



Another source of error is in the numerical differential of {(Ticr^j^-^) with re- 
spect to T. Since the derivative diverges at T^, a slight error in (cr^cr^_|_i) 
reduces the accuracy in Cv near T^. 

The numerical result shows that the DMRG method is quite efficient 
in the off-critical — high- and low-temperature — regions. On the other 
hand, it is difficult to maintain good numerical accuracy near the critical 
point, because long-range spin correlations prevent us from obtaining good 
renormalized transfer matrix, whose size (= 2m) is limited by computational 
restrictions. This convergence property of the DMRG method for the Ising 
model is in accordance with the fact that the DMRG method works more 
efficiently for ID quantum systems with a finite excitation gap.^^ 

I would like to thank S. R. White and M. Guerreo for helpful comments 
and discussions. I also thank Y. Akutsu, M. Kikuchi, K. Okunishi, and 
H. Otsuka for valuable suggestions. Most of the numerical calculations were 
done by using the super computer NEC SX-3/14R at the computer center of 
Osaka University. 
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Figure Captions 

Figs. 1. Renormalization processes, (a) We assume that the transfer matrix 
T(2^) can be renormahzed as the product form T^^^) in eq.(2). (b) The 
assumption is justified through the relation between T^^"* and j'j^^^^ 
in eq.(7). 

Fig. 2. NumericaUy calculated specific heat of the Ising model. The plotted data 
are obtained within 1024 iterations (2A'" = 2048) under the condition 
m = 60. The numerical error e in each data point is indicated by the 
symbols. 
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